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Abstract 

A hard sphere fluid confined by hard, structureless, and parallel walls is in- 
vestigated using a certain version of weighted density functional theory. The 
density profile, the excess coverage, the finite size contribution to the free en- 
ergy, the solvation force, and the total correlation function are determined as 
function of the slit width L for various bulk number densities pb- In quantita- 
tive agreement with rigorous results the present version of density functional 
theory yields a constant and large but finite number density profile for the 
limiting case that L is reduced to the diameter of the hard spheres. Within 
the Derjaguin approximation the results for the slit geometry allows us to 
obtain the solvation force between two large hard spheres immersed into a 
fluid of much smaller hard spheres. 
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I. INTRODUCTION 



The knowledge of the structural properties of fluids in confined geometries is important 
both for applied and basic research. For most applications one has to deal with an ensemble 
of interconnected pores with irregular sizes and geometries. This severely impedes a quanti- 
tative comparison between theoretical predictions and actual experimental data. Typically, 
in this case only general trends and spatially averaged quantities can be tested reliably. Con- 
sequently, in such systems many details of theoretical predictions for confined fluids remain 
unchecked. 

Therefore it is highly welcome that substantial experimental progress has been made 
to prepare model pores which consist of parallel plates whose surfaces are smooth both 
on atomic and mesoscopic scales and which are immersed into a fluid reservoir. Such well- 
defined systems are close realizations of corresponding theoretical models and can serve as 
a suitable testing ground for the behavior of the more complex systems mentioned above. 
Inter alia, by varying the distance L between the plates, one can study the crossover from a 
three-dimensional bulk system to a two-dimensional fluid. 

From an experimental point of view the structural properties of a fluid confined to this 
slit geometry can be probed on various levels. First, one can determine global properties 
such as the mean density in the slit and the excess density compared with a hypothetical 
bulk system of the same size. Second, ellipsometry and the reflectivity of light, X-rays or 
neutrons enable one to determine the density profiles normal to the slit surfaces. Third, 
atomic force microscopes allow one to monitor the solvation force acting on the two plates 
which reflects the change of free energy of the confined fluid as function of L. 

More resently, with the advent of powerful neutron and synchrotron sources a fourth 
component has been added to the spectrum of experimental techniques. The analysis of 
the diffuse scattering of X-rays and neutrons under grazing incidence gives access to the 
two-point correlation function of the confined fluid. Combined with the knowledge of the 
one-point correlation function, i.e. the density profile, this provides a deep insight into the 
structural changes of fluids induced by their confinement including lateral correlations (see, 
e.g., Refs. and references therein). 

The purpose of the present contribution is to provide a first step towards guiding such 
kind of experiments by calculating the two-point correlation function of a hard- sphere fluid 
between two hard walls (HSHW) based on a weighted cfensity-/unctional theory (DFT). Since 
this approach requires as a prerequisite the knowledge of the density profiles, we use this 
opportunity to compare our results for the one-point correlation function with those obtained 
previously by different techniques for HSHW; furthermore we pay particular attention to 
the limit of small values of L and to the discussion of the solvation forces. 

Our choice for this model system is motivated by its following virtues: 

(i) Due to its simplicity it is particularly well suited for comparisons with simulation data. 
Systems with soft repulsive or long-ranged attractive interaction potentials pose additional 
difficulties such as their perturbative treatment in analytic approaches and their unavoidable 
truncations in simulations. 

(ii) Within the framework of DFT long-ranged interactions between the fluid particles 
are typically incorporated by perturbation theory (see (i)), which needs as a prerequisite the 
results of the corresponding hard sphere reference systems. 
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(iii) The investigation of the HSHW model is not only an important step for the study 
of atomic fluids, but it is also an appropriate model for the description of other physical 
systems. Under favorable conditions certain colloidal particles between glas plates behave 
like hard spheres in a slit and can be investigated by means of video microscopy ||. Confined 
micelles represent another realization of this model system exhibiting substantial technical 
and biological interest 

(iv) The HSHW model is the simplest model which allows one to study a nontrivial 
dimensional crossover from d — 3 to d — 2. In this model the spheres loose one degree of 
freedom when the width L of the slit is reduced to the diameter of the hard spheres. This 
raises the question whether this system is purely two-dimensional and can be characterized 
as a hard-disc fluid, or whether the presence of the three-dimensional reservoir requires a 
different description. Since this is important for the interpretation of experiments with very 
narrow slits, in Sec. II we also introduce the density-functional theory of two-dimensional 
systems and provide a careful investigation of this limiting case in Sec. III. 

The HSHW model has already been investigated by a variety of techniques. For certain 
values of the chemical potential and of L simulations p|-|T3"l have provided density distri- 



butions and values for the solvation force. The same quantities and in addition the total 
correlation function have been studied in the framework of integral theories such as the 
Percus-Yevick approximation (PY) [p^-[T7|. Compared with these methods the DFT is com- 
putationally less demanding and also enables one to study the free energy of the system. 
Simulations face the difficulties that they are restricted to a few selected parameter values 
and that in the grand canonical ensemble fluids confined to narrow slits exhibit strong fluctu- 
ations |18[ . As far as the integral theories are concerned it is known that they do not capture 
interesting phase transitions such as wetting phenomena. Since for future work we are in- 
terested in them, we implement a specific form of DFT (weighted density approximation 
(WDA)) which is known to capture them. Thus in a later stage our approach will enable 
us to build on the present results for the description of the two-point correlation function 
close to such interfacial phase transitions. For these reasons it is worthwhile to analyze the 
HSHW model in terms of DFT. 

In Sec. [Tl] the DFT is introduced and the WDA used here is specified. In Sec. |TT| we 
investigate the limit L — > 2a. A thorough discussion of the density profiles (Sec. [IVD and of 
the correlation functions (Sec. |V|) follows. Section |VT] summarizes our main results. 

II. DENSITY FUNCTIONAL THEORY IN d SPATIAL DIMENSIONS 

A. One- and two-point correlation functions 

In thermal equilibrium the structural properties of a d-dimensional inhomogeneous fluid 
consisting of hard generalized spheres follow from the grand canonical partition function 



oo if 
ZiQ*, T; [V(R)]) = 1 + Y. J^m / ^ ■ ■ ■ ^ 

N=l ' J 



exp \B j d d R{[p, — V(R)]/5/v(R; {Rj})} — /3$({Rj}) 



(2.1) 
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as function of the chemical potential p and of the temperature T = l/iksP)] {Ri} = 
{Ri, . . . , Rat}. The particles are exposed to an external potential V^(R) which includes the 
confinement due to the container walls and thus limits the spatial integrations here and in 
the following. For the pair potential 

1 - 

*({R<}) = ? 5>(|Ri-R*|) (2.2) 
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with 



Z d describes hard spheres, discs, and rods of diameter a for d = 3,2, and 1, respectively. 
In terms of the number density operator /3jv(R; {Ri}) — X^i^Q^ ~~ ^i) the equilibrium 
density profile of the particles in the presence of the external potential V(R) is given by 

W (R)==^(R; {«,})> = -1 ^^™ . (2.4) 
The second derivative yields the total correlation function h(R±, R2) (pjv(R; {Ri}) = p(R)) : 

- - p = < p(Ri)p(R 2 ) > - < p(Ri) >< p(R 2 ) > 

=: h(R u R 2 )p d (R 1 )p d (R 2 ) + p d (R 1 )5(R 1 - R 2 ). (2.5) 

Within the framework of the density-functional theory the equilibrium density profile 
Pd(R) minimizes the grand potential functional [[HJ 

tt d ([p d (R)];T,p; [V(R)]) = F^(\p d (R)]; T) + F^([p d (R)\; T) 

- [ d d R(p-V(R))p d (R). (2.6) 

JR d 

The ideal gas contribution is given analytically by (A is the thermal de Broglie wave length) 
4%-(R)];r) = 4 / d*Rp d {R)\hi(p d Ql)A d )-l]. (2.7) 

The support of the trial function p~d(R) is that domain in M. d where the external potential 
V(R) differs from infinity; otherwise p d (R) = . The excess Helmholtz free-energy functional 
Fex 1 {[pd(R)]', T) is not known exactly and an appropriate approximation has to be chosen 
(see, c.f., Sec. II. B). Once the density profile p d (R) has been obtained by minimizing Eq. 
rrOD the direct correlation function 

cf (R..H.J MR)]) := ~0 ^l P f™ . (2.8) 

f)pd(Ri) op d (R 2 ) 
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and via the Ornstein-Zernicke equation 



^(R-i, R2) — ^(Rij R2) 

+ / rf d J R 3 4 2) (Ri,R 3 )p (i (R3)^(R3,R2) (2.9) 

the radial distribution function g d (R 1 ,R 2 ) = /id(Ri,R 2 ) + 1 are accessible. 

We shall compute all quantities for a slit consisting of two parallel structureless hard 
walls which are described by the external potential 

V(R) = \™ ' " <a ' Z > L -° . (2.10) 

Our choice of the origin and the width L are motivated by the comparison with an atomistic 
model of a slit. In this case the walls consist of two parallel semi- infinite crystals. The nuclei 
of the atoms forming the top layer of one of these crystals lie in a plane which is located at 
z = for the left wall and at z = L for the right wall. (We do not consider vicinal surfaces.) 
Between these two walls a fluid with a soft interaction with the substrate has a nonvanishing 
number density for < z < L. If the atoms forming the walls are replaced by hard spheres 
of diameter a and are smeared out in the lateral directions, a fluid with a hard sphere 
interaction is exposed to the potential defined in Eq. ( |2.10|) . Since this substrate potential 
is translationally invariant with respect to the lateral coordinates x and y (in d — 3), the 
density profile p 3 (R) depends only on the normal coordinate z as long as there is no freezing 



which leads to a periodic density variation also in the lateral directions ||20|| . The total and 
the direct correlation function depend on the normal distances z\ and Z2 from the wall and 
on the lateral distance r 12 = \r t — r 2 | where R = (t,z) = (x,y,z). (In the case of freezing 
the two-point correlation functions depend on ri 2 = ri — r 2 instead of |ri — r 2 |.) 



B. The linear weighted density approximation 

Since the exact expression for the Helmholtz free-energy functional is not known one 
has to resort to one of the approximations known in the literature |21| . Depending on the 



physical system and the quantities under consideration one chooses that approach which 
captures the essential features without being computationally too demanding. 

For a hard sphere fluid (d = 3) close to a single hard wall in a previous publication 
22|j we argued that the calculation of density profiles and of correlation functions can be 
carried out successfully using the linear weighted density approximation (LWDA) In 
this approach four weighted densities 

p„(Ri)=/ d 3 R 2 w u (\R 1 -R 2 \)p(R 2 ) (2.11) 
with normalized weights (00 = f cr 3 ) 
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1 , u = 

w u (R) = — B(a-R)t n 3 R\ 10q (2.12) 

V; I (1 + -)(1 - — ) , v = 1,2,3 
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are introduced. The excess free energy is a functional of these weighted densities, 

3 



+-\p(R)-p„(R)]f'M(R)) 



(2.13) 



where the functions 



— fvW 

V 



-16 + 4(1 - - log 1 -77 
V 

3(-16 + 26r7 - 7t7 2 ^ 
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2(1 -?7) 2 
- 68/7 + 25t7 2 



+ 3(l--)log(l-77) 



v 



V 



V 



-8(1 



n 



log(l - 77) 



1 

2 
3 



(2.14) 



depend on the dimensionless packing fraction 77 = uo p. By construction in the limit of a 



,(d=3) 



[P 



d=3 



and the 



homogeneous density distribution (p(R) = Pb) the LWDA free energy Fe X 

(2) 

corresponding direct correlation function c^ 3 (Ri,R2) reduce to the known Percus-Yevick 



(PY) bulk results ||24j| . This is important as we need a proper bulk limit in order to be able 
to describe correctly the influence of the walls. Furthermore this allows us to express the 
results of a slit, whose thermodynamic state is characterized by the intensive variable p, 
equivalently in terms of the density pb of a bulk fluid with the same chemical potential p. 
This facilitates the comparison with previous publications in which the results for similar 
geometries are expressed in terms of p\> [p2] ; in addition dependences on p& are easier to 
interpret than those on p. 

There seem to be only very few WDA which are specifically designed to describe an 
inhomogeneous hard disc fluid [ 25 1 . This dearth is tied to the fact that experimental results 
are rare, as it is difficult to realize a truly two-dimensional system experimentally and that 
the construction of many WDAs relies on the knowledge of an analytic expression for the 
bulk fluid free energy and of the direct correlation function. In three-dimensional systems 
of spheres the Percus-Yevick closure can be used, but there is no analytic solution thereof 
known for the two-dimensional case. (For an alternative approach see Ref. ||25|| .) Although it 
should be possible to construct a WDA also for this case following the concept of Curtin and 
Ashcroft ]26| which does not require analytic expressions for the bulk quantities, it is natural 
to analyze an alternative approach |27j which amounts to evaluate the excess free-energy 

(2.15) 

and which leads to the following approximate expression for the two-dimensional hard disc 
fluid in an area A: 

F ex\p2] = ^T J^dxdy !^p 2 (x,y)f' u (p u (x,y,z = 0)) 



functional Fet 3 ^[pd=3] (Eq. (|2.13|) ) of the hard sphere fluid for 

P3(R) = Pi(x,y) 8{z) 



i=0 
+ 



dz 



/,(p,(R)) - -p u (R)f' u (p u (R)) 



LWDA. 



(2.16) 
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The fact that this approximation originates from a theory designed for three spatial dimen- 
sions is especially apparent in the weighted densities 

p„(R)= / dx'dy'p 2 (x', y') w v {^\x - x'\ 2 + \y - y'\ 2 + z 2 ), (2.17) 
Jm 2 

which still depend on three coordinates. 



III. COMPARISON OF A HARD DISC FLUID WITH A HARD SPHERE FLUID 

IN A NARROW SLIT 



In order to assess the quality of the approximation leading to the free energy functional 
in Eq. ( |2.16| ) of a hard disc fluid we consider the special case of a homogeneous density 
distribution. By setting p 2 (x,y) = p 2 ,6 in Eq- (]2.15|) the excess free energy in Eq. ( [2.16| ) can 
be compared with the results of the scaled-particle theory (SPT) 



(3F^[p 2 , b ]=Ap 2 , b 



V2 



- ln(l - 772) 



SPT, 



(3.1) 



where A denotes the cross section of the slit and r] 2 = p 2b ja 2 the packing fraction. It turns 
out that for all densities the values of the LWDA free energy (Eq. ( |2.16| )) is higher than the 
SPT results and the difference increases with increasing density p 2>b . For p 2 ,b& 3 = 0.6 there 
is a deviation of about 17%. The pressure of the system is given by 



P 



dF^[p 2 



OA 



2 d(FW[p 2 ]/N) 
Pd~ 



dp 2 



(3.2) 



and can be compared with results of integral theories such as the hypernetted chain approx- 
imation (HNC) and simulations [2G]. The pressure calculated within the LWDA approxi- 
mation is comparable to that of the HNC results but higher than the one obtained from 
simulation data. Thus we conclude that the suggested functional in Eq. ( |2.16| ) is a reasonable 
but not very accurate approximation for a hard disc fluid. 

If the width of a slit filled with a hard sphere fluid is decreased one could be inclined 
to expect that in the limit L — > 2a (compare Eq. ( 2.1C )) the density profile reduces to a 
^-distribution as in Eq. ( |2.15| ) and consequently that in this limit the system is described by 
the density-functional theory of a two-dimensional system, e.g. by the one proposed above. 
However, if Eq. fl2.15p is inserted into the expression for the ideal gas contribution to the free 
energy functional (Eq. ( |2.7|) for d — 3), one obtains a mathematically ill-defined expression. 
Since this defect is not cured by the LWDA excess free-energy functional, a well-defined 
grand canonical functional Q[p] (Eq. (|2.6| )) can only be constructed, if p(R) remains a finite 
function even in the limit L —>■ 2a. Indeed Henderson [[30 
L := L — 2a the contact density is given exactly by 



has shown that in first order of 



p(z = a + ) = A 3 exp 



A 3 exp(/3/i) 




(3.3a) 



(3.3b) 
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This implies that in the grand canonical ensemble in the limit L — > 2a the fluid is squeezed 
out of the slit and that the number density < N > j A of the particles per area still contained 

inside the slit vanishes linearly as the width is decreased: < N > /A = a dzp(z) L —? 
p(z = a + )L. Due to this small number of particles per area the fluid behaves like an ideal gas 
and in zeroth order the density is determined by the Boltzmann distribution ( Eq. (|3.3b|) ). 
On the other hand the local number density is rather high because the value of the chemical 
potential p is imposed by the bulk reservoir. It is interesting to note that one obtains the 
same limit for rods confined to a finite line segment in the limit of a vanishing length of the 

t £-* i 



segment |31[| : < N > /L — > j-exp(pf3). 



In order to investigate the limit L — > within the framework of the LWDA the functionals 
in Eq. ( |2.13j ) and in Eq. (|2.7|) are simplified according to the following approximations. For 
small values of L the local density in the slit can be taken to be constant and equal to 
p(cr). Also the weights w v (\J\x\ — X2\ 2 + \y\ — 2/2 1 2 + \z\ — Z2\ 2 ) do not vary significantly for 
cr < z\ < L—a for v — 0, 1, 2, 3 and x\, X2, yi, y2, %2 £ K fixed. With the resulting functionals 



F^W)] = J2j A ^dy \^-p{a)Lfl{p v {x,y,z = 0)) 

(3.4) 

p„(R)= / dx'dy'p{a)Lw v {y/\x - x'\ 2 + \y - y'\ 2 + z 2 ), (3.5) 
Jr 2 



i=0 

+ / dz 



V,(p,(R))-^(R)/;(p,(R)) 



and 



(3F^[p(a)} =ALp(a) [ln(A 3 p(a)) - 1] (3.6) 
the grand canonical potential 

(3Q[p(a)} = (3 [Fg\p{a)\ + F^\p{a)}) - faALp(a) (3.7) 

is determined. Since p(a) remains finite, Lp(a) vanishes in the limit L — > and the functions 
/ in Eq. (p.4|) can be expanded into Taylor series yielding 



^ n[p(a)} = \a 2 ixp{a) 2 L 2 + p{a)L [ln(A 3 p(^)) - 1] - ^ L - 0. (3.8) 

Minimizing this expression with respect to the contact density p(a) leads to 

A 3 p(a) = exp(jJt/3 - p{a) na 2 L) (3.9) 



and via further expansion to Eq. ( p. 3a ). This renders the satisfying result that in the limit 



L — > up to first order in L the LWDA reduces to the exact result. 

Since as shown above the LWDA is capable of both describing reasonable well a hard 
disc fluid and reproducing correctly in first order of L the limit L — > 0, for a hard sphere 
fluid we are in the position to compare these two different physical systems within one and 
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the same approach. The excess free energy functionals in Eq. ( [2.161 ) and in Eq. ( |3.4| ) can 
be mapped onto each another using the replacement pib <-> Lp(a). This replacement only 
states, that the number densities per area have to be the same for both systems: 

/ dzp(z) = Lp(a) + 0(L 2 ) = p 2>b . (3.10) 

J a 

Obviously this mapping cannot be used to relate the ideal gas contributions of a hard disc 
fluid (Eq. ( p.7|) with d=2) and of a fluid in a slit (Eq. ( |3.6|) ). For them no simple map is 
found. Thus we conclude that the hard sphere fluid within narrow slits and connected to a 
reservoir does not resemble the genuinely two-dimensional hard disc fluid. If one wants to 
prepare a quasi two-dimensional system with a non-vanishing density p^ =2 one has to resort 
to the canonical ensemble, i.e., one has to restrict a fixed number of particles to a finite 
volume V = Aa. 

The mechanism, which leads to a finite density in the limit L — > 0, is revealed by Eq. fl3.8| ). 
The grand canonical potential is split into the excess part ^a 2 %p{a) 2 L 2 which captures the 
contribution of the interaction between the particles, the ideal gas part p(cx)L [ln(A 3 p(<j)) — 1] 
which mainly takes into account the entropy, and the part —p,Lp{a)f3 due to the chemical 
potential. As the excess part is quadratic in L it becomes less important as the width of the 
slit decreases and ultimately the ideal gas part determines the behavior of the system. Thus 
the entropy is responsible for the fact that the local density in the slit remains finite in limit 
L -> 0. 



IV. CHARACTERIZATION OF HARD SPHERE FLUIDS IN NARROW SLITS 

Since lateral ordering phenomena are beyond the scope of the present work, we focus 
on sufficiently small densities which are below the onset of such freezing transitions. In 
the homogeneous bulk freezing occurs at pbC 3 ~ 0.94 f32f , but already near a single wall 
prefreezing sets in at a slightly lower density |20]. Experiments in slits revealed a rich 



structure including phase transitions between different ordered states as the width of the slit 
varies. For these kind of systems canonical Monte Carlo simulation were able to reproduce 
the phase diagram satisfactorily ||35|| . By analyzing the total correlation function of Monte 
Carlo simulation data Chu et al. || have shown, that for L = 3a a hexagonal packing close 
to the walls occurs. Thus beside a complete freezing of the whole slit also lateral ordering 
in parts of the slit close to the walls seems possible. Therefore we have limited our present 
investigations to pb<J 3 = 0.68. For this density extensive Monte Carlo simulations [|TI|] gave 
no hints for an onset of lateral ordering. 



A. Properties of the density profile 

Using various mesh sizes (0.005<r . . . 0.05cr) for the integration the grand potential in 
Eq. ( |2.6|) was minimized for numerous slit widths (L = 2.001cr. . . 12a). If the width L is 
larger than 12<r, the resulting density profiles close to one of the walls agree well with those 



near the single wall of the corresponding semi-infinite system [22]. In this limit the wall 
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theorem (p(a + ) = f3P) is fullfilled. Figure |T] shows the density profile for two different slit 
widths {L = 5.1a and L = 3.8a). They are symmetric and exhibit a layered structure 
due to packing effects. Due to the absence of rigorous results the accuracy of our results 
can only be assessed by comparing them with simulation data. Surprisingly to our knowl- 
edge there is only a single publication, namely a molecular dynamics simulation (MD) |T3| . 
which allows us a quantitative comparison. In this case the chemical potential of the fluid 
is not known and one has to use a different parameter in order to be able to map the two 
approaches onto each another. Since the LWDA is not an exact theory, the degree of agree- 
ment depends on the choice of this parameter. Here we use the mean density p m (see, c.f., 
Eq. ( |4.3| )) of the particles in the slit, but we correct the value of the LWDA by a factor 
p(z = o~ + ) simulation / p(z = <j + )lwda = pr^- = 1-03 in order to take into account the fact that 
the LWDA slightly underestimates the pressure of the bulk fluid. Ppy denotes the pressure 



of a homogeneous LWDA fluid and Pes is the almost exact Carnahan-Starling pressure |36 



Figure 0(a) reveals a satisfactory agreement with the simulation data. Similar to the case of 



a hard sphere fluid close to a single wall of a semi- infinite system p2| , the first minimum is 
slightly too shallow, but the phases of the oscillations agree rather well. If one investigates 
the changes in the form of the density profile upon increasing the slit width L, one finds the 
following scenario: At very small widths the contact density p(a) is very high and the profile 
between the walls is almost flat. If L is increased the profile develops a U-shape with a single 
minimum at z = L/2 and both p(a) and p(L/2), which are shown in Fig. ^| as the full and 
the dotted line, respectively, decrease rapidly. As function of L the contact value attains a 
minimum at L — 2.39a and reaches again a maximum at L = 3.01er. As function of L the 
value of the density p(L/2) in the center reaches its first minimum at L = 3.07a and then 
increases slowly. At about L = 3.4a two local maxima in the density profile p{z) appear at 
approximately z = 2a and z = L — 2a (compare Fig. [TJ(b)) which are each approximately 
one hard sphere diameter a apart from the first layer at the distant wall. These two extrema 
merge into a single maximum when the slit width is further increased (L = 4.0a). It is diffi- 
cult to determine reliably the precise value of the slit width for which these two extrema can 
no longer be distinguished because it depends sensitively on the approximations entering 



LWDA p7\ . For L = 3.98cr the contact density attains its third maximum and the density 
profile has a pronounced W-like shape. For increasing L the peak in the center broadens 
and starting at about L = 4.5a it splits into two peaks located approximately at z = 2a 
and z = L — 2a. In between there is a local minimum, which deepens and is smallest for 
L = 5.1a (see Fig. |2[). The slit contains now four layers. If the width is further increased 
more layers are added by a similar mechanism. The extrema characterizing this process are 
given in Table |. The values of the contact density and of various other quantities which will 
be defined below (see Fig. |3|) oscillate as a function of L with a period of about a. Within 
each of these oscillations another layer is added to the slit. The density pl(L/2) at the center 
of the slit has a minimum if an even number of layers are in the slit and a maximum if there 
is an odd number. Therefore this quantity exhibits a periodicity of about 2a. 

A useful global description of the density distribution in the slit is given by the excess 
coverage 

T(L) := [ L dz[p(z)-p b ], (4.1) 
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and the related mean density 



Pm(L):=~ J dzp(z), (4.2) 

or 

p m (L)= Pb + T(L)/L. (4.3) 

As shown in Fig. ||] T(L) exhibits oscillations with a period of about a. Due to the factor 
1/L (see Eq. ( |4.3| )) the corresponding oscillations of p m (L) are less pronounced. 

In the limit L —>■ oo the coverage T(L) equals twice the coverage of a hard sphere fluid 
close to a single wall which has been discussed in Ref. [22]. For L larger than 5a the coverage 
differs only slightly from the limiting value r(oo): 

r(L ^7 \ {00) < 0.04, L > 5a. (4.4) 

r(oo) 

This specifies the range of validity for the well known approximation p m ~ pt + T(oo)/L 
I- ' 

In the limit L — > 2a the density p(z) is constant (Eq. ( |3.3b )) so that Eq. Q4.2p yields 

T(L - 2a)) = [L — 2a] [A -3 exp(/3/i) - p b ] + 0((^\ ) ^ 0. (4.5) 



B. Finite size contribution to the free energy and resulting solvation forces 

In the context of a slit geometry a particularly interesting quantity is the finite size 
contribution to the grand potential, 



[L - 2a] P PY ) 



(4.6) 



where Ppy is the pressure of a homogeneous bulk liquid at the same chemical potential p. 
In the limit L — > oo j(L) reduces to twice the surface tension of a hard sphere fluid close 
to the single hard wall of a semi-infinite system. Inserting Eq. ( |3.9|) into the expansions in 
Eqs. ( |3.S| ) and (|4.6| ) yields in the opposite limit L ^ 2a 



p(a) - P PY 



[L - 2a] + 0( 



L 



a 



(4.7) 



The behavior of ^(L) in the intermediate regime between these two limits is displayed in 
Fig. [| The difference j(L) — 7(00) decreases oscillatorily with increasing slit width. The 
maxima decay exponentially (~ exp(— 1.23L/cr)); their positions are given in Tab. Q. 

The force between the two plates is an experimentally accessible quantity 0. Based on 



thermodynamics this so-called solvation force per area f(L) is given as p9] 
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a-+oc A \dL 

Using Eq. ( f4.6| ) one obtains 



/(!) = - Urn 1 (^) -P t>) . (I.S) 

T,im,A 



so that 



j(L) - 7(00) = - [ L dL'f(L'). (4.10) 



Since 7(L -> 2a) = (Eq. (0)) it follows that 

POO 

7(00) = - / dL'f(L') = 2 lsl . (4.11) 



2(7 



Thus the surface tension 7^ = 7(oo)/2 of a hard sphere fluid close to a hard wall is a 
measurable quantity accessible to force measurements. Equation ( [4.1 1 ) remains valid even 



for softly interacting spheres close to a hard wall, but not if the hard wall is replaced by a 
soft substrate potential. 

The solvation force can also be expressed in terms of the difference between the contact 
density at the finite slit width L, p(a), and at infinite slit width, poo(cr), (see Ref. |I(J and 
the Appendix): 

f3f(L) = p(a)- Poo (a). (4.12) 

This difference is shown in Fig. [2| Minimizing the grand canonical potential within LWDA 
we obtain both the density profile as the minimizing function, which leads to the force via 
Eq. ([4.12| ) , and the value of the minimum Q [p] which yields the force using Eq. ( |4.8| ) . We find 
numerically that both routes lead to the same result. This can be anticipated because any 
WDA is thermodynamically self-consistent with respect to this relation (see the Appendix). 
Here it serves as a helpful check of the numerical calculations. 

For a physical understanding of the above results it is rewarding to consider the system 
depicted in Fig. |j. The solvation force per area f{L) is the net force exerted on the wall b 
and positive if it is directed outwards, i.e. to the right. For L > 2a it is shown in Fig. ^| and 
for a < L < 2a it is constant, 

f(L) = ~ Poo (a), a<L<2a } (4.13) 

because the particles on the right side of wall b exert the constant bulk pressure Ppy = 
ijPoc^o") to the left. Starting with the wall b at a position corresponding to a slit width L 
the work j(L) — 7(00) is gained (Eq. ( |4.10|) ) if one moves it to infinity (see Fig. ^j. If one 
now considers the particular case L = 2a~ the above process starts from a configuration 
involving only a single surface and leads to a configuration of three independent solid-fluid 
interfaces without changing the bulk contribution to the free energy of the total system but 
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increasing the surface contribution to the free energy by 2 / y s i. This provides a transparent 
interpretation of Eq. ( 4. lip because the integral over the solvation force is the work applied 



to the system during this process. In addition these considerations show that Eq. Q4.11D is 
valid in general and not only within LWDA. Together with the general relation in Eq. (|4.10|) 



this implies that the equation j(L — > 2a) = is also valid in general. Finally, since is 
negative |22] these arguments also tell that one gains work by moving the wall b from L = 2a 



to L = 00. Thus the generation of these two additional solid-fluid interfaces is favorable. 

The negative sign of 7^ may provoke the question whether the hard sphere fluid is 
actually stable against the spontaneous formation of cavities. The effect of the formation of 
such cavities on the density distribution of the fluid can be thought of to be same as the 
effect induced by the immersion of, e.g., a hard wall (Fig. |4j) with a minimum thickness a. 
The difference in free energy AF between the homogeneous and the corresponding perturbed 
system consists of the two surface contributions ^siA and the bulk free energy density times 
the excluded volume a A of the cavity: AF/A = 2j s i + Ppycr. Although 7^ is negative it 



turns out that AF is positive (41] so that the cavity formation is disfavored. 

The above discussion is concerned with the particular cases L = 2a and L = 00. For 
general L it is worthwhile to note that the extrema of j(L) correspond to the zeros of the 
solvation force F(L) which are documented in Table I. If in Fig. ^ the wall b is allowed to 
float freely, the minima of 'y(L) correspond to the most favorable slit widths. According to 
Fig. ^| the global minimum is located at L — 2.18a. Thus in thermal equilibrium the optimum 
configuration in Fig. f| corresponds to the case in which the walls a and b are separated such 
that a monolayer of hard spheres fits in between with a little bit space left. However, one 
should keep in mind that this statement is only valid if the mass of the piston is much 
larger than the masses of the hard spheres. Otherwise the fluctuations of the position of the 
piston must be treated on the same footing and together with those of the hard spheres. 
In this sense the above line of argument, i.e. that the equilibrium position of the piston is 
determined by the minimum of 7(£), corresponds to a Born-Oppenheimer approximation. 

At the width L = na just n — 1 spheres fit on top of each other into the slit. One 
may wonder whether this peculiar matching condition leaves a particular signature in the 
Independence of the various quantities studies above. The corresponding (d — 2)-dimensional 
problem of hard rods of length o confined to a segment of length L on a line can be solved 



exactly |!l|] and one finds in this case that the second derivative of the mean number of 
rods exhibits discontinuities at L — na whose magnitudes decrease for increasing values of 



L [51]. Since the additional spatial dimensions of a three-dimensional slit allow for an easier 
rearrangement of the spheres upon packing, we expect that in d = 3 these discontinuities 
are either smeared out or shifted to higher derivatives. Although in principle one should be 
prepared for the occurrence of such singularities in, e.g., ^(L) or f(L) they turn out to be so 
weak that they do not show up in our present LWDA approach on the scale of the numerical 
resolution we used. 



C. Derjaguin approximation for the force between large spheres 

So far our analysis has been confined to the study of the slit geometry which may be ap- 
plicable to force microscope measurements of confined colloidal particles whereby the solute 
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particles only contribute to the effective interaction between them; in the present context 
this effective interaction is approximated by a hard-core repulsion potential. However, it 
turns out that the results for the slit geometry can be even used to analyze this latter effec- 
tive interaction potential between large colloidal particles of radius R which are immersed 
into a solute composed of small particles with diameter a ||42|| . (For R —>■ oo this problem 
reduces to the standard slit geometry.) If the centers of the two large hard spheres are kept 
at a fixed distance 2R + h — a, in the limit R^> a the solvation force between them is given 



by the Derjaguin approximation [43 



POD 

f.{h)=nR dL'f(L') (4.14) 

Jh 

= ^[7(^-7(00)], 

where Eq. ( |4.11| ) has been used. The large spheres touch each other for h = a; in this case 
there is no small sphere between the two large spheres along the symmetry axis. According to 
Eq. ( ^4.14j ) the finite size contribution to the free energy of a slit of width h is proportional to 



the force between two hard spheres of radius R. For slit widths larger than 2a the free energy 
is given by Fig. |3] and for o < L < 2cr Eq. ( [4.13 ) is used in Eq. ( f4.14 ). The combination of 
these results leads to Fig. |[ 

The global minimum at h — a indicates that a strong depletion force will press the two 
spheres together if they touch each another. However, in order to find the thermodynamically 
most favorable separation one must consider the effective interaction potential 

W s = [ dtif s (ti) (4.15) 



between the two large spheres. This is shown in Fig. |5] as the dotted line. For low densities 
(per 3 < 0.2) this quantity has been investigated in the framework of an expansion into powers 
of p ||44|| . The present density-functional theory extends these results to higher densities. Inter 
alia, the interest in this effective potential arises from the question whether a binary mixture 
of hard spheres can exhibit flocculation. For such systems the PY theory rules out phase 
separation at all densities and size ratios [ 45| . However, modern integral theories indicate 



that phase separation can occur [46 



V. TWO-POINT CORRELATION FUNCTION 

The two-point correlation function of the HSHW model depends on the normal distances 
Zi and Zi from the left wall and on the lateral distance r 12 (see Subsec. II. A). The system 
is specified by the bulk density p^ of the corresponding homogeneous system with the same 
chemical potential and by the width of the slit L. Thus the correlation function depends 
on five independent variables. Within the context of a research paper a complete graphical 
account of the dependences on all five variables is not feasible. Therefore we have decided 
to discuss the general mechanism governing the behavior of the total correlation function 
on the basis of the Percus test-particle theorem (see the following paragraph) and to select 
the display of the dependence on r 12 (with z 1 = z 2 and fixed pb in all plots) for L fixed and 
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various values of z\ as well as for z\ fixed and various values of L. Here and in the following 
we discuss the case d = 3 only. 



In the present context the Percus test-particle theorem [|D| states that the product 
p(zi)g(r 12 , Zi, z 2 ) (see Eq. (|2.9|) ) for a hard sphere fluid in a slit of width L is equal to the 
one-point conditional density distribution p(Ri = (r 12 , £i)|R 2 = (0,z 2 )) of a hard sphere 
fluid exposed to an external potential consisting of a slit of width L and in addition of a hard 
sphere of diameter a whose center is fixed at R 2 = (0, z 2 ). For p^cr 3 = 0.546, L — 5.1, and 
z 2 = a this product is shown in Fig. [| In the limit r 12 —>■ oo it reduces to the density profile 
in Fig. [TJ(a) which corresponds to the same slit width and the same chemical potential. This 
comparison reveals that approximately p(Ri = (r 12 , £i)|R 2 = (0, z 2 )) is the superposition 
of the density profile of Fig. |l|(a) and the oscillatory density distribution around a fixed 
hard sphere placed in a previously homogeneous bulk fluid. Accordingly the coordinates Z\ 
of the maxima and minima of p(Ri|R 2 ) which are denoted in Fig. || by the dots and circles, 
respectively, almost coincide with those of the density profile of Fig. |I](a) and thus in Fig. 
|9] they line up nearly parallel to the wall. This general mechanism was also born out in 
previous analyses (PY-approximation for a fluid in a slit [[RJ; LWDA for a fluid close to a 
single wall p2[) and has proven to yield a roughly correct picture of the radial correlation 
function in confined geometries. 

It is rewarding to investigate for different slit widths L the radial distribution function 
9(^12, Zi, z 2 ) as function of r 12 , i.e. parallel to the wall with z\ = z 2 fixed. This reveals 
the influence of the distant wall on the correlation function close to the near wall. This 
dependence is of particular interest because it can be measured directly by using digital 
video microscopy [|]. Although such experimental data are not yet accurate enough to 
facilitate a quantitative comparison with theoretical results, we surmise that in the near 
future the experiments will be improved sufficiently. For a bulk density of p^er 3 = 0.683 the 
radial distribution function is shown for z\ = z 2 = a and z\ = z 2 = 1.5a in Figs. 0(a) and 
Fig. 0(b), respectively. Compared with the corresponding radial distribution function of a 
homogeneous bulk fluid with the same chemical potential (see the dashed-double-dotted 
lines in Fig. [?]) the amplitude of the oscillations is reduced for Z\ = z 2 = a but enhanced 
for Z\ = z 2 = 1.5a. Since already the increase of Z\ by the radius of a sphere alters the 
amplitudes so profoundly, we conclude that accurate measurements of parallel correlations 
require a spatial resolution in z-direction of about O.lcx or better. For increasing L the radial 



distribution function reduces rapidly to that of the corresponding semi-infinite system [22 



In Figs. 0(a) already for a slit width of L = 5a (L = 3a) for all values of r 12 the influence 
of the distant wall is no longer visible within the resolution of the plot. This is remarkable 
because for this width only 4 spheres fit side by side into the slit. 

In Sec. [TV] we have discussed the HSHW model for two sets of the parameters L and pb in 
terms of the one-point correlation function and related quantities. In Fig. [8] and |9| the same 
cases are now investigated in terms of the two-point correlation functions as function of r 12 
for various values of Z\ = z 2 and for p^a 3 = 0.546 and L/a = 5.1 in Fig. [8] (compare Fig. |l](a)) 
and for p&cx 3 = 0.683 and L/a = 3.8 in Fig. ^ (compare Fig. |T|(b)). The various values of z\ 
are chosen to coincide with the positions of the extrema in the corresponding density profile 
(see Fig. 0). For both systems in the case z\ = a the contact value g{ji 2 = a, z±, z 2 = zi) is 
strongly reduced as compared to the bulk value gpy(a), whereas for larger distances z\ these 
contact values are rather close to each other and to the bulk value. This is in accordance 
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with the results obtained for a semi-infinite hard sphere fluid near a single wall p2| |. A 



further analysis of these results reveals that for those values of Z\ which correspond to 
the minima of the density profile these amplitudes of the oscillations in ^(r 12 , Z\, z<i = Z\) 
are enhanced whereas for those values of z\ corresponding to the maxima of p(z) these 
amplitudes are reduced. This is surprising because the amplitudes of the oscillations in the 
bulk correlation function gpyix', Pb) decrease with decreasing pi,. Thus the natural attempt 
to approximate the radial distribution function such that in the case Z\ = Z2 it reduces to 
g(ri2, Zi, Zi = zi) ~ gpy( r = r \2\Pb = p( z i)) is unsuccessful because this approximation 
renders the opposite tendency of the actual oscillatory behavior at least for the values of 
the chemical potentials and slit widths considered here. On the basis of a known expression 
for the bulk correlation function g(r) there have been efforts ||47|j48| to construct a suitable 
ansatz for the radial distribution function for a hard sphere fluid close to a hard wall. Since 
such an ansatz does not incorporate the peculiar behavior of the actual radial distribution 
function described above, the reliability of these approximations is rather limited. 



VI. SUMMARY 

We have obtained the following main results for a fluid of hard spheres of diameter a 
confined to parallel and structureless hard walls at a distance L: 

(1) On the basis of the linear weighted density approximation (LWDA), which describes 
an inhomogeneous three-dimensional fluid in the grand canonical ensemble, we have con- 
structed an approximation for the two-dimensional analogon of a homogeneous hard disc 
fluid. This approximation compares favorable with simulation data. 



(2) As proven by Henderson [^0] in the limit that the width L of the slit is reduced such 
that it can accommodate at most a monolayer of the fluid {L — > 2a) the density profile 
approaches a large but finite constant value. This implies that in this limit the fluid is 
squeezed out of the slit. Up to first order in L = L — 2a we find that the LWDA reproduces 
this limit exactly. 

(3) The two-dimensional hard disc fluid and the three-dimensional hard sphere fluid 
confined to a narrow slit have been compared in the grand canonical ensemble. For finite 
chemical potentials the confined hard sphere fluid does not resemble the genuinely two- 
dimensional hard disc fluid. 

(4) The density profiles p(z) calculated within the LWDA compare satisfactorily with 
simulation data (Fig. [I]). The dependence of the contact density p(z = a) on the slit width 
L is close to that obtained from simulation data, although the amplitude of the oscillations 
in this dependence as obtained from LWDA is slightly smaller than that obtained from the 
simulations (Fig. H). 

(5) Both within the framework of exact thermodynamics and within LWDA the finite size 
contribution to the free energy j(L) (see Eq. ( f4.6| ) and Fig. |3[) represents the potential of the 
solvation force f(L). The minima of 'y(L) (see Table [I]) correspond to metastable distances 
between freely movable plates immersed into a fluid reservoir (Fig. Furthermore, within 
the Derjaguin approximation "f(L) renders the force between two large spheres suspended 
in a liquid of small spheres (see Eq. Q4.14D ). 
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(6) Within the LWDA and based on the above results we have determined the direct cor- 
relation function. By inverting the Ornstein-Zernicke equation the total correlation function 
has been calculated. We have discussed it within the framework of the Percus test-particle 
theorem (Fig. §). 

(7) The influence of the distant wall on the radial distribution function g(r 12, %\ = o~,z 2 = 
a) along the near wall decreases rapidly with increasing slit width L (see Fig. [5]). 

(8) For a fixed slit width and bulk reference density p& we have analyzed the dependence 
of the radial distribution function g(r 12 , Zi, z 2 = Z\) on z\. We find that for those values of Z\ 
for which the density profile p(z) exhibits a local minimum (maximum) the amplitude of the 
oscillations of this correlation function as function of r 12 is enhanced (reduced) compared to 
the corresponding bulk correlation function. 

A major advantage of the present density functional theory is that it is computation- 
ally much less demanding than integral equation theories or numerical simulations without 
loosing its competitiveness as far as the quantitative reliability is concerned. Furthermore 
density-functional theory yields relatively easy access to free energies and solvation forces. 
Therefore we are encouraged to extend this analysis to fluids governed by dispersion forces 
in order to refine the presently available ansatz for the two-point correlation functions in 
such systems P?|. 
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APPENDIX: 

Upon differentiating the equilibrium grand canonical potential in Eq. ( f2.6| ) one obtains 
by using the chain rule 

sa _ .,„, . f *u \ SFlp] ,,. wm ,„\ 6p(R') } 



SV(R) ry ' J M d Up(R') J 5V(R) 

For the equilibrium density distribution the expression within the curly bracket vanishes. 
This is true even for approximate expressions for the functional F[p] = Fid[p] + -F e x[p] like, 
e.g., the one used for the LWDA. According to Eq. ( j2.10|) the external potential V(R) 



depends parametically on L so that with Eq. (|Alj) one has 



(m\ = r j sn 8v(r) 



dexp(-/3V(R)) 
dL 



/ d d Rp(R)(~eMPV(R)) 
~ [ al d Rp(R) eMPV(R))-^-(Q(z - a) + 6(L - a - z) - I) 
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-A±{p{z)eMPV{z))) z=L _ a 



-A±p(L-a) = -A±p(a). (A2) 



The last but one step in Eq. ( |AJ ) is based on the fact that the product p(z) exp([3V(z)) is 
continuous as function of z so that the value of this product at z = L — a can be obtained 
by considering the limit z — > L — a — where V(z) = 0. This renders the equivalence between 
the two definitions in Eqs. fl4.8Q and ( |4.12j) both for exact and for approximate expressions 
for F[p]. 
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FIGURES 



FIG. 1. Density profile of a fluid of hard spheres with diameter a between two hard walls 
located at z = and z = L (compare Eq. ( p.lO| )) according to the LWDA density functional theory 
(full curve) for the slit widths L = 5.1a (i.e. L = 3.1a) (a) and L = 3.8a (i.e. L = 1.8a) (b). In 
(a) the corresponding bulk density is chosen as pi,a 3 = 0.546 to allow comparison with molecular 
dynamics simulation data In (b) the bulk density is p^a 3 = 0.683 as in all remaining figures. 



FIG. 2. Within the LWDA the contact density p(a) (full line) and the density p{L/2) in the 
center of the slit (dotted line) of a hard sphere fluid are shown for a bulk density of p^a 3 = 0.683. 
The slit width L varies between 2.1a and 8a. In the limit L — > oo the contact density approaches 
the constant value Pooi^a 3 = 3.84. In the limit L — > 2a the two densities p(a) and p(L/2) approach 
the same limit A -3 exp(p/kT) as given by Eq. ( 3.3b| ). The dots denote results of grand canonical 
Monte Carlo simulations Ol. 



FIG. 3. The coverage T(L) and the finite size contribution to the free energy j(L) (Eq. ( |1.6[ )) 
of a hard sphere fluid between two flat hard walls for various slit widths L = 0.1a . . . 8a at a bulk 
density of p^a 3 = 0.683. In the limit L — > oo the coverage approaches r(oo) = — 0.92cr 2 and j(L) 
reaches the value (3^(oo)a 2 = —1.95. 



FIG. 4. Schematic plot of a model system consisting of a planar hard wall (a) on the left side 
and a hard piston (6) which is in front of and parallel to the left wall. It can be fixed at different 
distances L. The hard sphere fluid between a and b and to the right of b are connected to the same 
grand canonical reservoir and thus are in equilibrium with each another. The force per area and 
UbT acting on the piston is given by /3f(L) = p(a) — Poo{&) and is plotted in Fig. [2| as the full line. 



FIG. 5. The force f s (h) between two fixed hard spheres of a large radius R immersed into 
a fluid of hard spheres of diameter a <C R as function of the minimum distance h between the 
surfaces of the large spheres (see Eq. Q4.14D ). The centers of the large spheres are at a distance 
2R + h — a. The dotted curve represents the effective interaction potential between the two large 
spheres. 
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FIG. 6. Within LWDA the conditional singlet density p(Ri[R,2) = p(z%) <?(ri2, zi, z% = a) of 
a hard sphere fluid in a slit of width L = 5.1a at a reference bulk density of pb<r s = 0.546 in the 
presence of a hard sphere of the same diameter whose center is fixed at R2 = (r%2 = 0, z% = cr). 
(This position is marked by a cross.) Due to g(ryi —* 00, Z\, Z2) = 1 one recovers for ryi — > 00 Fig. 

a) which corresponds to the same values of L and pb- The dots and circles in the contour line 
plot at the bottom of the figure denote the positions of the local maxima and minima, respectively. 
The dashed line indicates the excluded volume due to the test particle; it is not an isodensity line. 
The value of the singlet density varies by an amount of O.lcr -3 between neighboring contour lines. 
The contour lines are shown only for values less then 1.0er~ 3 . (In the actual calculations a mesh 
size of 0.02o" in the z- and the redirection has been used.) As can be seen from the contour lines 
the perturbation of the density distribution in the slit due to the presence of the fixed sphere at 
R-2 = (0, a) dies out for ri2 <L 3<r or z\ > 3a. 

FIG. 7. The radial distribution function 5(^12,-21,^2) of a hard sphere fluid for various slit 
widths L at the reference bulk density p^a 3 = 0.683. The lateral distance ri2 varies between 1.2a 
and 3.0a for z\ = Z2 fixed with z\ = a in (a) and z\ = 1.5a in (b). In the limit r%2 — > a the 
correlation functions increase rapidly to values 2.3 . . . 2.4 in (a) and 2.8 ... 3.1 in (b) depending on 
L. Note that we have used the same scales of the axes in (a) and (b) in order to facilitate a direct 
comparison of the two cases. In both plots the dashed-double-dotted curves denote the correlation 
function of the corresponding homogeneous bulk fluid within the PY-approximation. In (a) for all 
values of ri2 the correlation function cannot be distinguished from its semi-infinite form p2| for 



FIG. 8. The radial distribution function g(ri2, Z\, Z2 = z{) of a hard sphere fluid in a slit of 
width L = 5.1a at a bulk density p^cr 3 = 0.546. The decay of the correlations parallel to the walls 
is shown for various values of z\ = Z2- According to Fig. 1(a) z\ = l.Ou corresponds to the contact 
with the left wall, Z\ = 1.64cr to the first minimum, z\ = 2.2a to the first maximum, and Z\ = 2.55<7 
to the midpoint of the density profile of this system. 

FIG. 9. The radial distribution function g(r±2, z\, Z2 = z\) of a hard sphere fluid in a slit of 
width L = 3.8cr at a bulk density pf,a 3 = 0.683. As in Fig. || the decay of the correlations parallel 
to the walls, i.e. as function of ri2, is shown for various values of Z\ = Z2- According to Fig. 1(b) 
Z\ = 1.0a corresponds to the contact with the left wall, z\ = 1.46cr to the first minimum, z\ = 1.78<r 
to the first maximum, and z\ = 1.9a to the midpoint of the density profile of this system. 
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TABLES 



TABLE I. Characteristics of a hard sphere fluid between two hard walls of width L at a bulk 
density p;,<7 3 = 0.683. The table lists those slit widths for which the contact density p(cr), the 
density p(L/2) at the center, the coverage L(L) (Eq. (|4.l| )), and the finite size contribution to the 
free energy 'y(L) (Eq. (|4.6D) attain their extrema. 



width L/a leading to minima of 

p{a) 2.39 3.41 4.42 5.43 6.43 7.42 

p(L/2) 3.07 5.12 7.12 

L(L) 2.63 3.62 4.60 5.60 6.60 7.60 

7(L) 2.18 3.12 4.22 5.21 6.21 7.20 

width L/a leading to maxima of 

p(a) 3.01 3.98 4.95 5.95 6.93 

p(L/2) 4.12 6.15 

2.06 3.10 4.10 5.10 6.09 7.08 

7(L) 2.76 3.74 4.72 5.71 6.71 
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